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When  the  flow  over  a  submerged  ,  round,’ upright  cylinder,  situated  in  a 


large  ocean,  is  forced  by  a  train  of  plane  waves,  linear  theory  (Yamaiuuro, 
1981)  shows  that  the  response  can  be  abnormally  large  for  certain  forcing 
frequencies.  The  aim  of  this  paper  is  to  present  a  weakly  nonlinear  theory, 
where  wave  interactions,  arising  from  the  quadratic  terms  in  the  free-surface 
boundary  conditions,  can  yield  abnormally  large  responses. 

A  specific  interaction  will  be  considered  between  a  flow  at  a  subharmonic 
frequency  and  a  flow  at  the  driving  frequency.  The  reason  for  considering 
such  an  interaction  derived  from  a  consideration  of  some  experimental  results 


of  Barnard,  Pritchard  and  Provis  (1981) 
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SIGNIFICANCE  AND  EXPLANATION 


Observations  of  unusually  large  response  in  the  wave  records  at  Macquarie 
Island,  in  the  ocean  south  of  Australia,  have  been  explained  as  being  due  to 
the  excitation  of  nearly  resonant  trapped  modes*  This  phenomenon  has  moti¬ 
vated  much  experimental  investigation  of  trapped  modes  around  a  submerged  up¬ 
right  cylinder  in  water  of  constant  depth.  Theories,  based  on  the  small- 
amplitude  approximation,  lead  to  predictions  of  resonances*  However,  some 
experimental  work  indicated  the  presence  of  subharmonic  components  in  the 
wavefield.  This  paper  develops  a  possible  explanation  for  this  by  showing 
that  nonlinear  interactions  could  yield  resonance  of  subharmonic  modes. 
Moreover,  a  particular  set  of  conditions  is  presented,  in  which  the  amplitudes 
of  plane  waves  incident  on  the  obstacle  are  significantly  magnified  above  it 
by  subharmonic  resonance* 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report* 


WEAKLY  NONLINEAR  INTERACTIONS  AND  WAVE-TRAPPING 


Yuriko  Renardy 


$1.  Introduction. 

If  a  train  of  plane  waves  of  a  certain  frequency  is  incident  on  a 
submerged  round  sill,  situated  in  a  large  ocean  (Figure  1),  linear  theory 
(Yamamuro,  1981)  predicts  that  the  overall  amplitudes  over  the  sill  may  become 
much  larger  than  those  of  the  deeper  ocean.  Such  a  phenomenon  may  be  called 
"wave-trapping"  or  "near-resonance",  referring  to  the  unusually  large 
response.  The  purpose  of  this  paper  is  to  examine  a  contribution  of  nonlinear 
effects  to  the  wave-trapping  phenomenon. 

Let  the  fluid  over  the  submerged  sill  be  denoted  by  and  the  fluid 
outside  this  region  by  D2 .  The  total  domain  for  the  flow  will  be  denoted 
by  D  «  D(  u  D^.  The  boundary  conditions  at  the  free  surface  are  nonlinear. 
However,  the  amplitudes  of  the  motion  in  D2  are  assumed  to  be  small  enough  so 
that  those  conditions  may  be  linearized.  An  investigation  is  made  of  a  wave¬ 
trapping  phenomenon  in  which  'small-amplitude'  waves  enter  D.,  and  become 
magnified,  so  that  in  this  region  some  nonlinear  terms  are  included  in  the 
free-surface  conditions.  A  possible  mechanism  of  magnification  is  a  wave 
interaction  arising  from  the  quadratic  terms  in  the  free-surface  boundary 
conditions  (Phillips,  $3.8),  and  it  seems  appropriate  to  restrict  attention  to 
these  in  an  initial  investigation.  Therefore,  if  0  is  the  driving 
frequency,  attention  will  be  limited  to  interactions,  in  D^ ,  with  waves  of 
frequencies  0/2  and  20. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041.  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS-7927062. 


Figure  1 


The  structure  of  the  paper  is  as  follows. 

§2,  }3.  In  §2,  the  equations  governing  the  flow  are  presented.  An  attempt  to 
satisfy  them  by  a  modal  decomposition  of  the  velocity  potential  in  the  radial 
variable  leads  to  an  infinite  matrix  equation.  The  organization  of  the  solu¬ 
tion  method  for  this  set  of  equations  is  described  in  $3.  A  linear  calcu¬ 
lation  ( Yamamuro,  1981)  shows  that,  in  order  to  make  the  velocity  continuous 
across  the  sill-edge,  the  wave-field  necessarily  contains  modes  that  decay 
away  from  the  sill-edge  and  are  named  "decaying  Diodes".  Resonance  of  such 
modes  will  not  be  considered  here.  Hence  the  solutions  exhibiting  large 
response  are  assumed  to  be  combinations  of  eigenstates,  each  of  which  consists 
of  a  principal  part,  which  is  wavelike,  together  with  associated  decaying 
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ill.  In  order  to  calculate  the  amplitudes  of  the  wave  motion,  it  is 
necessary  to  solve  a  'homogeneous'  problem  in  D  ,  namely  a  linearized  problem 
with  no  incident  waves.  This  is  investigated  in  §4.  The  problem  posed  in  D2 
is  solved  from  an  existing  linear  theory.  The  flow  structure  in  D2  is  then 
used  to  generate  a  boundary  condition  at  the  sill-edge  r  ■  1  so  that  the 
problem  is  reduced  to  a  boundary- value  problem  in  D^.  It  is  found  from  a 
separation  of  variables  that  an  eigenstate  in  Dj ,  periodic  in  time  with  a 
complex-valued  frequency  ft  ,  consists  of  a  'wavelike'  mode  and  an  infinite 
number  of  'decaying'  modes.  In  addition,  there  are  an  Infinite  number  of 
complex-valued  coefficients  to  be  determined  from  the  boundary  conditions. 

The  eigenfrequencies  ft  and  the  coefficients  have  been  determined 
simultaneously  by  finding  the  zeros  of  an  expression  in  the  depth  variable.  A 
collocation  method  is  used  and  consists  of  applying  that  expression  at  N 
values  of  the  depth  variable  and  including  the  first  N  modes  of  each  eigen¬ 
state.  This  yields  an  N  *  N  matrix,  whose  zeros  can  be  found  by  the 
following  method.  The  condition  number  of  the  matrix  is  computed  over  a  grid 
of  complex  frequencies  and  those  with  relatively  large  condition  numbers  are 
taken  to  be  approximations  to  the  ft's.  This  method  was,  however,  time- 
consuming  and  an  alternative,  non-standard,  method  was  devised.  The  method  is 
an  iterative  scheme,  based  on  the  smallness  of  the  response  of  the  decaying 
modes  compared  with  that  of  the  wavelike  modes. 

i ii  The  solution  to  the  homogeneous  problem  provides  a  set  of  orthogonal 
functions,  a  combination  of  which  can  be  used  to  form  an  expression  for  the 
surface  elevation.  The  total  sill  solution  is  then  written  as  a  suitable 
combination  of  the  orthogonal  functions,  superposed  on  the  linear  solution. 

In  )5,  an  example  of  such  a  nonlinear  interaction  is  constructed.  The  example 
chosen  here  displays  subharmonic  resonance  and  was  motivated  by  previous  work 
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on  edge  waves  (Guza  &  Davis,  1974;  Minzoni  &  Whitham,  1977}  Rockliff,  1978)* 
The  interaction  involves  three  modes:  two  cos  20  modes  at  the  forcing 
frequency,  and  a  cos 9  mode  at  half  that  frequency.  It  is  found  that  near¬ 
resonance  occurs  for  two  ranges  of  forcing  frequencies.  One  range  occurs  near 
an  eigenfrequency  of  the  cos20  mode,  and  the  other  occurs  near  twice  an 
eigenfrequency  of  the  cos 9  mode. 


I-  • '  ■.  vr  •»  =  ■*-»  %r 


§2.  Formulation  of  the  problem. 

It  is  assumed  that  the  motion  is  inviscid  and  irrotational •  The  domain 

of  the  flow  above  the  sill  (r  <  1)  is  denoted  by  and  that  outside  the 

sill  region  (r  >  1)  by  Dj*  The  amplitudes  of  the  flow  in  D2  are  assumed 

to  be  sufficiently  small  so  that  the  free-surface  boundary  conditions  may  be 

linearized  and  applied  at  an  equilibrium  level.  Such  a  linearization  is  not 

assumed  in  The  geometry  of  the  problem  is  shown  in  Figure  1. 

Cylindrical  coordinates  r,9,z  are  used:  r  is  measured  outward  in  units  of 

the  sill  radius  a  ,  and  z  is  measured  vertically  upward  in  units  of  the 

undisturbed  water  depth  d  above  the  sill.  The  depth  outside  the  sill  is 

denoted  by  D  «■  d/6  .  The  velocity  potential  $(r,9,z,t),  which  is  in  units 
2 

of  d  cr,  satisfies  Laplace's  equation,  namely, 

2 

$  +  —  $  +  — —  $  ■  0  .  (2.1) 

rr  r  r  2  SB  2  zz 
r  d 

Let  n(r,9,t)  denote  the  surface  displacement  and  x  denote  the  horizontal 
coordinate  in  the  direction  9*0.  A  train  of  plane  waves  is  assumed  to  be 
incident  from  the  positive  x-axis  and  is  represented  by  the  real  part  of 


_  -i(kx  +  at) 

n  =  In^  e 

where  k  is  the  positive  real  root  of 

2 

(KD/a)  tanh  (KD/a)  ■>  Dd  /g. 

The  free-surface  boundary  conditions  are: 

at  z  *  h(r,9,t),  0  <  9  <  2*,  r  >  0,  a  kinematic  condition 

♦  -  n/0  +  (d2/a2)(V  *)  •  (V  n) 
z  n  n 

where  V  =  (3/3x),  3/3y),  and  a  dynamical  condition 

n 

d2°(<^tt  +  (*/d,V  +  H* )  +  V2(u**  V(u*»  u*))  -  0 


(2.2) 


(2.3) 


where 


u  -  d2o($  /a,  t9/ar,  ♦  /d). 
—  r  v  z 


(2.4) 


(2.5) 


(2.6) 


•.  SW’  W>".  *•' 
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The  small-amplitude  approximation  of  condition  (2.4)  to  second  order  is  found 


(Phillips#  $3.1)  to  be#  for  z  -  0#  0  <  9  <  2t#  r  >  0# 

♦tt  +  (g/d)*z  ■  -h  <*tt  +  ( g/d ) *z )  -  -^(u**  u*)/d2o.  (2.7) 

The  possibility  of  a  resonance#  in  which  the  right-hand-side  becomes  of  the 
same  order  as  the  combined  expression  on  the  left-hand-side#  will  be  inves¬ 
tigated  for  the  sill  region.  The  possibility  of  a  similar  kind  of  resonance 
through  condition  (2.5)  will  not  be  considered.  Under  these  assumptions#  the 
free-surface  boundary  conditions  to  be  used  are:  for  z  «  0# 

T-11  "Sa  f^tt  +  *  21  J/d2®  for  D1 

+  (g/d)*  -  <  (2.8) 

“  *  for  D2 


and 


n  ■  0*  for  D  .  (2.9) 

t  z 

As  mentioned  in  $1,  only  flows  of  frequencies  0,  2c  and  c/2  will  be 
investigated,  so  that  condition  (2.8)  can  be  expressed  in  the  following  form: 


*tt  +  9/d  ♦* 


f , (r#8)  e"i0t+f2(r,9) 


for 


-  V2  iot 


+f,(r,8)e 


-2iot 


+  *  for  D 


1 


(2.10) 


Here#  the  asterisk  denotes  the  complex  conjugate  of  preceding  terms.  The  f^ 
(i  “  1 #2,3)  are  complex  functions  of  *  and  are  determined.  Once  *  is 
known  and  substituted  into  equations  (2.8)  and  (2.9). 

The  boundary  condition  on  rigid  surfaces  is  one  of  zero  normal  velocity# 

i .  e . 

*  ■  0  for  z  -  -1#  r  <  1  and  z  -  -1/5,  r  >  1  ,  (2.11) 

z 

and 

♦r  -  0  for  r  -  1,  -1/5  <  z  <  -1.  (2.12) 
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At  the  sill-edge,  the  velocity  is  assumed  to  be  continuous 
that  is  equivalent  to  the  $  and  being  continuous  there, 
distances  from  the  sill,  the  wavefield  is  assumed  to  consist  of 
wave,  and  waves  that  either  decay  or  radiate  outward. 
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a  condition 
At  large 
the  incident 
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§3.  Form  of  'resonant'  solutions 


The  problem  formulated  in  §2,  with  (2.8)  and  (2.9)  as  the  free-surface 

conditions,  is  essentially  a  superposition  of  two  problems.  Hence  the 

solutions  are  expressed  as  ♦  *  $_  +  .  The  <p  satisfies  the  linearized 

Ii  N  Ju 

free-surface  conditions,  namely, 

4  +  (g/d)4>  =0  at  z  =  0  ,  0  <  9  <  2*,  0  <  r  <  “  ,  (3.1) 

L  tt  Lz 

and  is  forced  by  the  plane  waves  represented  by  (2.2)  and  (2.3).  The  ♦  is 

constructed  to  make  the  total  velocity  potential  $  satisfy  the  nonlinear 

free-surface  conditions  (2.8)  and  (2.9),  and  is  not  forced  by  the  plane 

waves.  Both  *  and  ♦  satisfy  the  conditions  (2.11)  and  (2.12)  at  the 
L  N 

solid  walls.  The  $  interacts  with  ♦_  through  the  free-surface  conditions 

N  L 

and  the  solutions  of  interest  are  those  in  which  ,  becomes  comparable  with 

M 

*L* 

The  linear  solution  <b  has  been  calculated  (Yamamuro,  1981)  by  a 

Ii 

separation  of  variables  for  the  regions  0  <  r-  <  1  and  r  >  1,  and  the 

velocity  is  made  continuous  throughout  the  flow. 

The  <frN  are  expressed  as  a  linear  combination  of  orthogonal  functions 

that  satisfy  equation  (2.1)  and  conditions  (2.8),  (2.11),  (2.12),  and  the 

radiation  condition  that  the  flow  be  non-growing  at  large  r  .  Let  the  part 

of  $N  in  be  denoted  by  and  that  in  D2  by  $N2»  $N2  is  to 

satisfy  linearized  conditions  and  to  consist  of  flows  of  frequencies  0,  c/2 

and  20.  Separation  of  variables  then  yields 

00 

♦  _  =  e  I  cosm9{B  ^ ( aX  r/D)coshXA(6z+1 ) 

N2  _  mo  mu  u 

m=0 

00 

+  I  B  K  (aX  r/D)cosX  (6z+1)/K  (aX  /D)}  +  *  (3.2) 

_  mn  m  n  n  ran 

n=«1 

+  (similar  expressions  for  flows  at  frequencies  0/2  and  20} 


where  (±Xq,  ±1^#  iiX^,***}  are  the  roots  of  the  dispersion  relation  (Davis 
&  Hood,  1976) 


X  tanh  X  =  do  /g. 
.(1) 


(3.3) 


Notations  for  the  Bessel  functions  (x)  and  K^x)  are  those  used  by 

Abramowitz  and  Stegun  (1972).  The  complex  coefficients  are  yet  to  be 

determined.  However  a  boundary  condition  for  is  first  constructed  by 

Nl 

the  elimination  of  the  Bm  . 

mn 

At  r  =  1,  the  condition  that  $  and  3$  /3r  be  continuous  is: 

N  N 


3<£  /3r  = 

n2/ 


3  $  /3r 

NV 


for  -  1  <  z  <  0 
for  -1/6  <  z  <  -1, 


(3.4) 


and 

4>  =  <{>  for  -1  <  z  <  0  (3.5) 

N2  Nl 

The  expression  (3.2)  is  substituted  in  equation  (3.4).  In  what  follows,  the 

discussion  is  focused  on  flow  at  one  of  the  three  frequencies  (o,  a/2  or 

20)  denoted  by  <*>.  The  orthogonality  of  the  set  {cosh  X^(6z+1), 

cos  X  (6z+1):  n  =  1,2,***}  over  -1/6  <  z  <  0  is  used  and  integration 
n 

over  z  yields  the  Bmn  in  terms  of  3  <^>N ^ / 3 zr  at  r  =  1.  These  equations 
are  then  used  in  equation  (3.5)  to  eliminate  the  Bmn»  A  boundary  condition 

_iU)t 

results  for  $  (z),  the  coefficient  of  cos  m0e  in  <$>  at  r  =  1: 

m  Nl 


3$ 

$  (z)  =  f°  (z*  )k  (z,z'  )dz'  for  -1  <  z  <  0  ,  (3.6) 

m  ^  or  w 

where 

H  (aX/D)  cosh  X  (6z+1)  cosh  X  (6z'+1) 

v  iz  \  =  _ _ _ _ _ _ 5 _ 9 _ 

V  ,z  H  ’ (aX/D) (aX/D)  h(X  ) 

m  0  0  o 

(3.7) 

®  K  (aX  /D)  cos  X  (5z+1)  cos  X  (6z'+1) 
r  m  n _ n _ n _ _____ 

+  L.  K'(aX  /D)(aX  /d)  h ( i X  ) 
n=  1  m  n  n  n 
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Thus,  the  remaining  problem  is  to  solve  for  using  equations  (2.1), 

(2.11)  and  (3.6).  The  details  are  given  in  §4.  For  each  of  the  frequencies 

w  ( a ,  a/2  or  20),  the  corresponding  flow  in  is  constructed  as 

follows.  The  spatial  variation  of  ♦  is  constructed  to  be  that  of  the 

-ifl^t 

solutions  4  (r,z)cos  m0e  +  *  which  satisfy  the  "homogeneous"  problem 

mn 

for  ♦  ,  i.e.  condition  (2.8)  is  replaced  by  the  linearized  form.  The  set 

N  l 

(r,0);  m  fixed,  n  ■  1,2,***)  can  be  shown  to  be  orthogonal  for 
mn 

0  <  r  <  1  by  the  application  of  Green's  theorem  (Jeffries  &  Jeffries, 

§5.081)  to  the  domain  D..  Thus  ♦  can  be  written  as 

*  N  i 


♦hi  *  \  008  m0  t.{amn*'iat/2  *amn(r'2) 

m»0  n*1 


-10  t 

where  ♦. _ (r,z)cos  m©  e  ,  (r,z)cos  m&e  and 

amn  Pmn 


(3.8) 


-i«  t 

*  (r,z)cos  m©  e  are  the  solutions  to  the  homogeneous  problem  for 

iron 

<*>  equal  to  o/2,  o  and  20  respectively.  The  complex  coefficients  a  , 

mn 

6  and  Y  are  determined  by  condition  (2.8),  in  which  the  orthogonality 
mn  mn 

property  of  those  functions  is  used.  For  example,  the  «  are  determined  by 

mn 


a 

mn 


(-o2/ 4  +  n2  )  JV 

'  amn  '  amn 


(r,0)rdr 


«  JU 

0 


(r,0) 


/  2*  f  (r,0)cos  m0  d0 
0  2 _ 

/21Tcos2m0  d© 


rdr. 


L  0  J 

Similar  equations  hold  for  B  and  Y  .  These  equations  in  which 

ran  ran 

and  interact  will  be  referred  to  as  ’interaction’  equations. 


(3.9) 


* 
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It  can  be  seen  that  the  'response'  lo^j  is  inversely  proportional  to 

\-d*/4  +  fl2  I  so  that  in  actual  computations  involving  a  particular 
1  ctmn 

geometry,  only  the  one  a^,  whose  is  the  closest  to  0/2,  is 

expected  to  be  important.  Similarly,  at  most  one  8^  and  one  are 

expected  to  be  significant  enough  for  inclusion  in  the  interaction 
calculations. 

We  note  in  passing  that  the  Longuet-Higgins  eigenfunctions,  referred  to 
as  'free  modes'  in  his  paper  (1967),  cannot  be  used  here  for  the  construction 
of  <frN1  because  then  the  range  of  integration  on  the  left  hand  side  of 
(3.9)  would  extend  to  infinity  and  the  integrand  would  grow  almost 
exponentially  with  r  . 
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at  z  -  0,  ♦  +  (g/d)$  -  0 

tt  z 

at  z  -  -1,  ♦  ■  0 

z 


§4.  The  homogeneous  problem  on  D^« 

This  problem  consists  of  equation  (2.1) ,  with  the  boundary  conditions; 

(4.1) 

(4.2) 

at  r  *  1,  condition  (3.6)  holds  and  w  is  assumed  to  be  the  frequency  of 
the  flow  in  D2< 

Separation  of  variables  yields  solutions  of  the  form 
0(r,z)co8  m6  e  where 


(4.3) 


4(r,z)  ■  A  J  (akr/d)cosh  k(z+1) 
m  zn 


+  £  A  1  (ak  r/d)cosk  (z+1)/I  (ak  /d) 

.  mn  in  n  n  in  n 

n=1 


and  Am,  \in*  k,  kn  jj  are  to  b®  determined.  Equations  (4.1)  and  (4.2)  yield 
a  dispersion  relation  between  k,  kft  and  fl  ,  namely 


r*. 

{  }  tanh  {  }  *  dfl  /g. 


(4.4) 


ik 


ik 


The  A m,  \m'  *  and  are  relate<*  through  condition  (3.6),  which  yields 

the  following  equation; 


A  L(k,z)  +  I  A  L  (k  ,z)  ■  0  for  -1  <  z  <  0 


m 


n*1 


mn  n  n 


(4.5) 


where  the  operators  L(k,z)  and  Xin(kn,z)  are  defined  in  the  appendix.  The 
wave  numbers  X  and  X^f  appearing  in  the  operators  L  and  Ln,  correspond 
to  the  fixed  frequency  w  of  the  flow  in  t>2»  They  are  known  through  the 
dispersion  relation 


{  }  tanh  {  }  *  Dw2/g 


iX 


iX 


(4.6) 


n 
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where  u  »  o,  0/2  or  20 


A  matrix  equation  that  expresses  the  in 

terms  of  An  and  also  yields  the  A  ,  is  constructed  by  satisfying  equation 
(4.5)  at  N  values  of  z  and  neglecting  the  coefficients  A^  for  n  >  N. 
Convergence  with  N  was  checked  numerically.  One  way  of  obtaining  the  0's 
is  to  search  through  the  complex  plane  by  computing  the  condition  number  of 
the  matrix  over  a  grid  in  the  complex  plane.  If  a  grid  point  yields  a  near¬ 
singular  matrix/  then  that  complex  number  can  be  taken  to  approximate  an  A  . 
An  alternative,  less  time-consuming  method  is  now  described. 


§4.1.  Iterative  scheme. 

The  scheme  is  constructed  to  take  advantage  of  the  largeness  of  lAjJ, 

representing  the  response  of  the  wavelike  modes,  as  compared  to  lA^I 

representing  the  response  of  the  decaying  modes,  and  of  the  property  that 

lA^I  is  smaller  for  larger  n  j  i.e.  the  wavefield  contains  little  of  the 

modes  that  decay  very  fast  away  from  the  sill-edge.  Equation  (4.5)  is 

multiplied  by  each  of  the  functions  in  the  set  {cosh  k(z+1),  cos  k  (z+1): 

n 

n  =  1,2,***}  and  integrated  over  -1  <  z  <  0.  Since  the  elements  of  that 
set  are  orthogonal  to  each  other,  the  resulting  equations  take  the  forms 
for  m  =  0,1,2,*** 


A  x(k,A,A  )  +  l  A  X  (k  ,A,A  )  *  0 
m  n  ms  s  s  n 

s*1 


and  for  p  =*  1,2,  *  *  * 


Y(k  ,k  , A, A  )  +  l  A  Y(k  ,k  ,A,A  )  -  -a  Y(k,k  ,A,A  ). 


mp  p  p 


8*1 


ms  s  p 


m 


(4.7) 

(4.8) 


The  functions  X,  X^  and  Y  are  defined  in  the  appendix. 


The  iteration  proceeds  as  follows*  The  first  iterate  for  k,  denoted 
by  k^,  is  calculated  from  the  following  reduced  form  of  equation  (4.7),  in 
which  all  decaying  modes  are  neglectedt 

J  (ak/d)H*(aA/D)f(k)«h(A)A  -  k  J'(ak/d)H  (aX/D) (g(A,k)]2  -  0  (4.9) 

m  m  in  in 

where  X  is  known.  The  corresponding  eigenfrequency  is  calculated  from 

k(0)tanh  k(0)  -  dfl(0)2/g,  (4.10) 

after  which  the  k^°^'s  are  calculated  from  equation  (4.4).  Next,  the 
equations  (4.8)  are  used  to  express  A^/J^  for  n  ■  1,2, •••  .  Then,  on 
using  these  relations  to  eliminate  the  A^'s,  equation  (4.7)  takes  the  formt 

x  (k)  +  SXh(X)x  (k,k  )H'(aX/D)  -  0  (4.11) 

1  2  n  m 

where  x^(k)  represents  the  left  hand  side  of  equation  (4.9)  and  contains  no 
decaying  modes,  x2(k,kn)  involves  the  decaying  modes,  and  the  notation  is 
defined  in  the  appendix. 

The  nth  iterate  k^  (n  ■  1,2,***)  is  calculated  by  a  Newton's 
method  from  equation  (4.11)  in  which  the  term  x2(k,kn)  is  calculated  at  the 
known  (n-l)th  iterate,  i.e., 

x,(k(n))  +  «  Ah(A)  x  (l^H"*1  ^  .k^11*"1  ^  )H'  (aX/D)  -  0  .  (4.12) 

The  corresponding  ft*n^  and  k^n)  are  then  calculated,  as  for  the  zeroth- 
iterate  from  equation  (4.4). 

A  numerical  check  of  the  scheme  was  performed  as  follows.  The  ft's  were 
compared  with  the  eigenfrequencies  for  the  domain  D  denoted  by  ftQ  because 
they  were  expected  to  have  similar  values.  The  ftp's  generalize  the 
frequencies  of  the  'free  modes'  investigated  by  Longuet-Higgins  (1967)  using 
shallow-water  theory,  and  were  calculated  in  a  similar  way  to  the  ft's.  The 
only  difference  was  that  in  equation  (4.11),  X  was  also  an  unknown.  There¬ 
fore,  at  each  step  of  the  iteration,  the  two  equations,  (4.11)  and 

k  tanh  k  -  5  X  tanh  X  ■  o  ,  (4.13) 

were  solved  simultaneously  for  k  and  X  by  a  Newton's  method. 


•I 

i 


§5.  Example  of  a  near-resonance . 

A  particular  set  of  conditions  in  which  the  foregoing  theory  yields  near- 
resonance  will  be  presented.  In  order  to  simplify  computations#  the 
parameters  5  and  d/a  will  be  chosen  to  be  small.  The  smallness  of  d/a 
ensures  the  smallness  of  the  effect  of  the  decaying  modes  in  the  flow  in  D<j 
so  that  in  the  'interaction'  equations#  such  as  (3.9)#  the  decaying  modes  will 
be  assumed  to  be  negligible.  However#  the  decaying  modes  will  not  be 
neglected  in  the  computation  of  the  ft's  since  these  are  required  to  a  high 
order  of  accuracy.  Furthermore#  the  smallness  of  6  ensures  that  some  of  the 
ft's  will  have  very  small  imaginary  parts#  so  that  if  the  flow  is  forced  near 
such  an  eigenfrequency,  near-resonance  is  possible. 

The  experimental  scales  of  Barnard#  Pritchard  and  Provis  (1981)  were 
examined  for  the  presence  of  near-resonant  nonlinear  interactions  but  were 
found  not  to  yield  them.  However#  a  choice  of  scales  which  do  ares 
d  -  2^,  d/a  -  .005,  «  «  .002,  | n 1 1  -  .01. 

In  this  case,  the  maximum  amplitudes  of  the  decaying  modes  is  at  least  an 
order  of  magnitude  less  than  that  of  the  wavelike  mode.  Although  these  scales 
are  unusual#  there  may  be  other  realistic  combinations  where  the  nonlinear 
theory  is  applicable. 

An  interaction  of  three  modes  will  be  considered:  the 

(cosde  iot^2  +  *)  an<j  (Cos  20e  i0t  +  *)  modes  in  and  the 

N 

(cos  2@e  +  *)  mode  in  i  .  In  the  sill  region#  these  are  represented  as 

L 


follows: 


♦  *  A  J  (akr/d)  cos20  cosh  k(z+1)e  +  *  +  decaying  modes  (5.2) 

Ii  2  2 


$  ■  a  cosB  e  (aVr/d)cosh  v(s+1)  ♦  *  ♦  decaying  modes 

N  1 

♦  B  cos  20  e”*0tJ2(aJJr/d)cosh  M(z+1)  +  *  +  decaying  modes 


(5.3) 


where 


•J 

l 


k  tanh  k  *  do2/g,  v  tanh  v  -  d  ft2/g,  (5.4) 

M  tanh  M  ■  dfl2/g. 

(1^  represents  the  8  defined  in  $4  which  lies  closest  to  0/2  for  the 
cos2$  mode  and  ^  Is  the  fl  which  is  closest  to  0  for  the  cos20  mode. 
The  free-surface  boundary  conditions  (2.8)  and  (2.9)  yield 

o  R  -  i  o*  (A  V  +  BV2)  (5.5) 

BR2  -  ia2v3  (5.6) 

where 

R  *  (d [02/g  -  d02/4g)cosh  v  /^J2(avr/d)rdr  (5.7) 

11  0  1 

R  ■  (dfl2/g  -  do2/g)cosh  y  J1  J2  (a»ir/d)rdr.  (5.8) 

2  2  o  2 


The  functions  V1 ,  V2  and  V3  are  defined  in  the  appendix.  The  response  for 
the  linear  forcing  A2  can  be  calculated  from  a  method  described  in  (Yamamuro 
1981).  A  trivial  solution  is  a  =  0  and  B  *  0.  The  questions  to  be 
resolved  are  whether  there  are  any  other  solutions,  and  if  so,  under  what 
conditions. 

Eliminating  B  from  equations  (5.5)  and  (5.6)  yields 

a/a*  *  i  (A2V1/V2  +  i  a2  V^/JR^))  (5.9) 

Let 


le1*. 


10. 


-V2V3/(RlV 


If  le 


and 


Then 


i02 

9  *  i  a2V2/Ri  “Isle 
(cos  t  (cos21  -  1  +  |g|2)  ^2]/|f| 


(5.10) 

(5.11) 


and 
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(5.12) 


♦  “  ^2  ( ®2  +  sin”1  ( |  <*| 2  |f  Isin  ^/Igl)). 


Next,  B  is  evaluated  via 


B  -  i  a2V3/R2.  (5.13) 

In  order  that  there  be  nontrivial  solutions,  two  conditions  must  be 
2 

satisfied.  First,  |a|  must  be  positive  and  from  (5.11),  the  conditions 
that  must  be  met  are: 

(a)  cos2S1  -  1  +  |g|2  >  0  , 

(b)  if  cosQ^  +  (cos29^  -  1  +  |g|2)  0  , 

there  is  at  least  one  nontrivial  |a|.  if 

cos^  -  (cos291  -  1  +  |g|2)  ^  >  0  , 

then  there  are  two  solutions  for  |a|.  Secondly,  equation  (5.12)  shows  that 
2 

I  a |  | f | / 1 g|  must  be  less  than  or  equal  to  1  . 

If  is  close  enough  to  <J,  then  k  is  approximately  so  that  the 

total  velocity  potential  in  the  sill  region  is,  approximately, 

♦  ■  a  cos  9  e  (aVr/d)cosh  v(z+1) 


+  (B+A  )cos20e  (akr/d)cosh  k(z+1)  +  *. 

d  dm 


(5.14) 


This  approximation  may  be  used  for  the  present  example.  In  this  case, 
computations  revealed  two  ranges  of  forcing  frequencies  0  ,  in  which  near¬ 
resonance  occurs.  One  range  lies  near  2^  and  the  other  is  near  0^  *  In 
most  of  these  ranges,  the  wave  amplitudes  were  calculated  to  be  rather  high  so 
that  instabilities  may  occur,  after  which  the  present  theory  may  not  be 
applicable.  However,  at  the  upper  end  of  the  range  near  2ft^,  the  amplitudes 
were  found  to  be  small  enough  so  that  the  present  steady-state  theory  might  be 
observable  in  practice. 
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S K  ■».-v  -  .  i  :=_- . 


u 


1 

'  i 


r* 

4 


■'Mr,  4-  .-  • 


L  (k  ,z)  ■  cos  k  ( z+1 )  - 
n  n  n 


H^(aX/D)  cosh  X( Sz+l)g(X,ik  ) 
m  n 


/ 1  >  • 

„  H  (aX/D)  X  h(X) 
L  m 


"  K  (aX  /D)  cos  X  (6z+l)g(iX  ,ik  ) 

+  l  —2 - £ - , - El. _ £ _ E_ 

K'(aX  /D) X h{iX  ) 


p*i  op  p 


k  I'(ak  /d) 
n  ra  n 

6  I  (ak  /d) 
m  n 


k  J* (ak/d) 
o 

H-(aX/D) (g(X,k)]2 

in 

5 

H'(aX/D)  X  h(X) 
o 

•  K  (aX  /D)[g<iX  ,k)]' 

+  y  »  p  p 

\  KMaX“/D)  X  h(iX  ) 
P“*  m  p  p  p 


-I' (ak  /d)k 

X  (k  ,X,X  )  -  -  2  ■  E-—  £ 
p  p  n  I  (ak  /d)  6 
m  p 


H  (aX/D)g(X,ik  )g(X,ik  ) 
ra _  p  p 

H * ( a VD )  X  h(X) 


;  v*V0,9(lV,tl9<1ViV 
*  ll  K;(*VDI  X„  h(1V 


Y(k, k  i X , X  ) 
P  n 


k  J'(ak/d) 
n 


H  (aX/D)g(Xfik  )g(X,k) 

m  p 

H' (aX/D)  Xh(X) 


ra 


-  Km(«>n/D)9(l>n,lkD),(lXn,10 

,c;<aVD>  x»h<iv 


I* (ak  /d)k 

Y(k  fk  ,X,X  )  «  f(ik  )  -  -?-  £■  ■  <6 
p  p  n  p  1  (ak  /d)5 


m  p 


Hm(aX/D)[g(X,ik  )]' 
m  p 

H' (aX/D)  Xh( X) 


ra 


•  K  <*\./D)  lg( iX  ,ik  )]' 

V  S  n  n  p 

K'(aX  /D  X  h(iX  ) 
n- 1  ra  n  n  n 
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k  I'(ak  /d)  H  (aX/D)g(X,ik  )g(X,ik  ) 
>.i.  *  \  \  t  sms  w  p  s 

Y(VV  'V - y  r(ST7S)  H*  (aX/D)  Xh'(X) 

n  3  m 


•  K  (aX  /D)g(iX  ,ik  )g(iX  ,ik  ) 

T  P  t  t  P  Z  8 

A  K'(aX./D)  x  hCiX  ) 

t"  1  mu  t  u 


x  (k,k  ) 
2  n 


kJ  (ak/d)  •  K  (aX  /D)[g(iX  ,k)]‘ 
n  r  p  p  p 

S  K* (aX/D)  X  h(iX  ) 

P“1  m  p  p  p 


*  A  I'(ak  /d)  k  H  (aX/D)  g{X,k)g{X,ik  ) 
T  jojp  p  jd _  p  p _  _  p 

^  A  I  (ak  /d)  H"  (aX/D)  X  h(X) 

p*1  m  m  p  I  p 


-  g(iX  ,k)g(iX  ,ik  )  X  (aX  /D) 
Y  n  n  P  p  n 

L4  X  h(iX  )  K'(aX  /D) 
n-1  n  n  m  n 


_  2  _ 

r  *  -  [-2  v  ainh  v(-aZk  sinh  k  +  kZ  cosh  k)  +  k  sinh  k(-  —j  v  sinh 


♦  v  8inh  v)]  x  j^Ur 


+  a2(  %  k  cosh  k  v  cosh  v  J!  (^■^)<J,(^x^)jl(~r^)rdr 
4  .la  la  2  a 


♦  t  *  sht? 


+  V2k  sinh  k  v  sinh  v  / 

0 


where  k  tanh  k 
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^2  identical  to  but  with  M  instead  of  k  ,  where  y  tanh  y 
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